Jahn- Teller distortions and phase separation in doped manganites 



A. O. Sboychakov, K.I. Kugel, and A. L. Rakhmanov 
Institute for Theoretical and Applied Electrodynamics, 
Russian Academy of Sciences, Izhorskaya Str. 13/19, Moscow, 1254-12 Russia 
(Dated: 6th February 2008) 

A "minimal model" of the Kondo-lattice type is used to describe a competition between the 
localization and metallicity in doped manganites and related magnetic oxides with Jahn- Teller ions. 
It is shown that the number of itinerant charge carriers can be significantly lower than that implied 
by the doping level x. A strong tendency to the phase separation is demonstrated for a wide 
range of intermediate doping concentrations vanishing at low and high doping. The phase diagram 
of the model in the x — T plane is constructed. At low temperatures, the system is in a state 
with a long-range magnetic order: antiferromagnetic (AF), ferromagnetic (FM), or AF-FM phase 
separated (PS) state. At high temperatures, there can exist two types of the paramagnetic (PM) 
state with zero and nonzero density of the itinerant electrons. In the intermediate temperature 
range, the phase diagram includes different kinds of the PS states: AF-FM, FM-PM, and PM 
with different content of itinerant electrons. The applied magnetic field changes the phase diagram 
favoring the FM ordering. It is shown that the variation of temperature or magnetic field can induce 
the metal-insulator transition in a certain range of doping levels. 
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I. INTRODUCTION 

The effect of electron correlations on the properties 
of different materials is currently among the most burn- 
ing problems of the condensed matter physics. As a 
rule, strong electron correlations are accompanied by 
the formation of nanoscale inhomogeneous statesi. Such 
inhomogeneities have been already studied for several 
decades. In particular, they were widely discussed for 
high- Tc superconductors^ and heavy fermion materials^. 
In the recent years, they attract a special attention ow- 
ing to the discovery of the colossal magnetoresistance ef- 
fect in manganites (the nature of which is believed to 
be closely related to inhomogeneous structures^). The 
inhomogeneities manifest themselves in other magnetic 
materials such as cobaltitesi^, nickelate9& and also in low- 
dimensional magnets^. All these systems are character- 
ized by a strong interplay of spin, charge, and orbital 
degrees of freedom leading to rather rich phase diagrams. 

One of the first most spectacular examples of such 
kind of inhomogeneities is the formation of ferromagnetic 
(FM) droplets (magnetic polarons or ferrous) in antifer- 
romagnetic (AF) semiconductors at low doping levels as 
well as FM spin polarons in the paramagnetic state^iS. 
These examples correspond to the case of so called elec- 
tron phase separation caused by self-trapping of charge 
carriers, which change their local environment. In addi- 
tion to such a small-scale phase separation, in mangan- 
ites, as well as in other compounds exhibiting first order 
transitions (e.g., between FM and AF phases), there also 
arises the phase separation of another type related to 
rather wide region where different phases coexist. An 
example of such large-scale phase separation is the for- 
mation of rather large FM droplets with the size of the 
order of (100-1000) A inside the AF matrisiMi. At 



higher doping levels close to the half-filling, there ap- 
pears one more threshold for the phase separation in the 
system corresponding again to the formation of ferro- 
magnetic droplets, but now in a charge-ordered insulat- 
ing matrisi^. The interaction of spin, charge, and or- 
bital degrees of freedom can also lead to the formation 
of stripe structures instead of droplets at high content of 
the alkaline-earth elementi^. In manganites, owing to the 
strong electron-lattice interactions, such structures are 
related to the lattice distortions and can be observed by 
the electron diffraction and the low-angle neutron scat- 
teringi^. 

Both analytical and numerical studies in various mod- 
els related to the strongly correlated electrons exhibit a 
pronounced tendency toward phase separation in a wide 
range of temperatures and electron or hole concentra- 
tions. Among these, we can mention s — di^, t — 
Hubbardii, and Falicov-KimbalU^ models. 

The theoretical models usually imply that the num- 
ber of charge carriers introduced by doping is equal to 
the number of itinerant electrons which take part in the 
formation of nanoscale inhomogeneities. However, the 
comparison of experimental data with the theoretical re- 
sults suggests that such an approach is insufficient^^ and 
the number of self-trapped carriers can significantly differ 
from the doping levelSSiSi. 

Here, we analyze the model proposed in our paper 
Ref. 123. which relates the doping level and the number 
of charge carriers. We take into account the Jahn- Teller 
(JT) nature of magnetic ions, which could give rise to the 
localization of charge carriers at the lattice distortions. 
We introduce this localization effect to the Kondo-lattice 
model in the double exchange limit. Beginning from sem- 
inal paper, Ref 2^ the role of JT distortions was widely 
discussed in literatur o^i^^i^^ . In particular, such distor- 
tions were taken into consideration in Refs. in the 
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analysis of the phase diagram of doped manganites. How- 
ever, none of these papers dealt with the phase separa- 
tion phenomena. In our paper, the main emphasis is 
made on the studying the phase diagram in the doping- 
temperature (x — T) plane and its evolution under effect 
of the applied magnetic field. 



II. THE MODEL 



First, let us note that the doped manganites, 
(Re^t^Mn^j:^)(A2+Mn;^+)03", are the compounds with 
mixed valence. Here Re is a trivalent rare-earth element 
and A is a bivalent alkaline-earth element. The ions 
Mn'^+ and Mn'*+ have 3d* and electron configura- 
tions, respectively. In such compounds, manganese ions 
are located in the centers of Og octahedra. In a reg- 
ular octahedron, a five-fold degenerate 3c? level is split 
into triple- and double-degenerate levels, t2g {dxy, dy^, 
dzx) and Cg (dj.2_y2, dy,z2_r^), respectively. The t2g level 
lies lower than eg level. Manganese ions are character- 
ized by the strong Hund's rule coupling giving rise to 
the parallel alignment of intra-atomic electron spins. So, 
the spins of t2g electrons form local spin S — 3/2. In 
a regular MnOg octahedron, Mn'^+ ion has one electron 
at double-degenerate eg level22i. According to the Jahn- 
Teller theorem, the latter configuration is unstable and 
the degeneracy is lifted by a distortion of the octahedron. 
The Mn^+ is not a Jahn- Teller ion and Mn'^+Oe octahe- 
dron remains undistorted. The distortion of Mn'^+Og oc- 
tahedron leads to the energy lowering by cjt- In doped 
manganites, the eg electron can hop from Mn'^'^ to Mn^+ 
ion producing a gain in the kinetic energy due to electron 
delocalization. Therefore, an electron can lower its en- 
ergy either due to the Jahn- Teller-induced localization 
at distorted octahedra or by the delocalization related to 
the inter-atomic hopping. The strong Hund's rule cou- 
pling favors the hopping of an electron when its spin is 
parallel to the spin of core {t2g) electrons. This is the 
origin of the well-known double-exchange mechanism of 
ferromagnetic interaction between the localized spins^S. 

Therefore, it is natural to assume that the eg electron 
can be either localized due to Jahn- Teller distortions with 
the energy gain ejT {I electron) or to decrease its energy 
due to band broadening {b electron)2SiSL^. So, there ex- 
ists a competition between the localization and delocal- 
ization. Such a system with localized and band electrons 
can be analyzed using the Hubbard Hamiltonian taking 
into account the electron-lattice interaction, the Hund's 
rule coupling, and the exchange interaction between core 



electronsSS 
H 



He\ + Haf + HjT + Hcl-cl , 

^ ^ ^ ^ ] i^nmO-liacr'^mba + h.C. 
(nm) aba 

'^naa ('^S„)^^, anaa 



(1) 



(nm) 

3n('^ ^ab]^nba 



n aba 



K 



l-cl 



n 

Ui \ ^ U2 \ ^ 

/ ^ ^naa^naa i / ^ '^naa'^ 



In this Hamiltonian, aj^^g. , Unaa are creation and anni- 
hilation operators for eg electrons at site n with orbital 
index a (82^ — or — y^) and spin projection a, Sn is 
a local spin of t2g electrons. Below we will consider Sn as 
classical vectors, cr are the Pauli matrices and Q2in Qsn 
are normal modes of vibration of MnOg octahedron. The 
symbol (nm) denotes the summation over nearest sites. 
The electron part -ffei of Hamiltonian |^ includes the ki- 
netic energy of eg electrons and the Hund's rule coupling 
between the spins of eg and t2g electrons. Haf is the 
antiferromagnetic (J' > 0) exchange interaction between 
local spins. The iJjT term takes into account interac- 
tions between eg electrons and vibrational modes for the 
MnOg octahedra, here K is the elastic energy and g is the 
electron-lattice coupling constant. The on-site Coulomb 
repulsion Hd-d includes the terms corresponding to eg 
electrons at the same and different orbitals, where the 
bar above a or a means not a or not cr, respectively. 

We consider the limit Jh ^ 00 characteristic of man- 
ganites. In this case, the spin of eg electrons is par- 
allel to Sn and we can eliminate the spin indices by 
the transformation of a^aa to operators Cna with spin 
projection -1-1/2 onto the direction of Sn accompanied 
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by the transformation of hopping amplitude; 
*nm cos(i^nm/2), whcrc COS f nm = SnSm/S"^. In addi- 
tion, we assume that t^ln — > for / electrons, which pro- 
duce maximum splitting of eg level with the energy gain 
—g'^/2K, whereas b electrons with non-zero hopping inte- 
grals tnm produce smaller distortions of MnOg octahedra. 
Preliminary calculations for the case inm ^ demon- 
strated that the results are not significantly affected if 
the 'localized' band is much narrower than itinerant one. 
Therefore, in this paper, we consider the limiting case of 
zero hopping integral for 'localized' electrons. Hamilto- 
nian then reads 
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where ribn = cJ^Cn and — l^ln are the numbers of 
b and / electrons at site n, cjj, c„ and Zjj, are the 
creation and annihilation operators for the b and I elec- 
trons, respectively, and is the chemical potential. The 
first three terms in HcS correspond, respectively, to the 
kinetic energy of b electrons, JT energy of localized elec- 
trons, and on-site Coulomb repulsion between b and I 
electrons. The last term in Hos is the Heisenberg anti- 
ferromagnetic exchange between local spins. The effec- 
tive on-site Coulomb repulsion U in Eq. |(2Jl can differ 
from U2, but has the same order of magnitude (~ 5eV). 
ejT ~ g^f^K is the JT energy gain for I electrons count- 
ing from the center of b electron band. The number of 
localized, n/, and band, Tif,, electrons per lattice site obeys 
an obvious relation nh + ni = 1 — a;, where x is the doping 
level. 

In our paper, we limit ourselves by the case of the 
large on-site Coulomb repulsion U, which strongly sup- 
presses double occupancy of the site. Moreover, at large 
U, the characteristic time of existence for a 'double elec- 
tron state' is of the order of l/U <C where $7 is the 
characteristic JT phonon frequency {h — 1), which is of 
the order of Debye frequency. Therefore, the adiabatic 
approximation for the Jahn- Teller distortions is applica- 
ble. As a result, the JT term in Eq. Q is determined 
only by localized electrons. 

Hamiltonian ^ was analyzed at zero temperature in 
Ref. The homogenous ferromagnetic and antiferro- 
magnetic states as well as the phase-separated FM-AF 
state were studied. The effective parameters t and ejx 
were considered to be independent of the densities of 
band rib and localized ni electrons. In the present analy- 
sis, we also neglect the dependence of parameters ejx and 
t on rifc and n;, but take into account the temperature 
dependence of hopping integral t, which can be rather 
strong due to polaron band narrowing22iS. Following 
Refs. iSasl we write the expression for the hopping inte- 
gral in the form 

t(r)=ioexp(^--^^75^j . (3) 

In this expression, A is the dimensionless electron-phonon 
couphng constant. From Eq. it is clear that the hop- 
ping integral t{T) decreases with temperature. Even if 
we take into account the finite bandwidth for 'localized' 
electrons, it will also decrease with T. However, the pa- 
rameters describing this behavior could be different for 



different bands. We cannot assert that at high tempera- 
tures the ratio of widths for the narrow and wide bands 
would be larger or smaller than at low temperatures. In 
this paper, we consider the temperature range, which is 
relatively small as compared to characteristic value of Q. 
Thus, t(T) does not vary by orders of magnitude and the 
bandwidth ratio remains small. 

To study the effects of temperature, we use the mean 
field (MF) approximation. For this purpose, we make a 
decoupling procedure in the first term of Hes, Eq. ||2J), 
in the following way. The values of c|,Cm and i?nm = 
\/(S'2 -t- S„Sm)/252 = V(l + enem)/2 Can be written 
as 

cl^C^ = (4cm) + S (c|^Cm) , "i^nm = (l?nm) + <5l9nm , 

where angular brackets mean thermal averaging, Sn = 
Se^, and e„ is the unit vector. Omitting the products 
proportional to 6 (c|jCm) S'&nm, we write the first term in 
HcS as 

-t 4cm- ("^nCm) ' 

(nm) (nm) \ / 

where t — t{T) (1 + e„em)/2y Note that in homoge- 
nous state, (i?nm) does not depend on indices n and m 
(for sites m nearest to the site n). Now, the effective 
Hamiltonian can be represented as a sum of electronic 
and magnetic parts 

n 

Hc\ = X! '^n'^'" ~ '^JT X! + ^ X! "'n^lbn, (4) 
(nm) n n 

= - ^ [(4c^) (M„™ ~t)- J'S^B^e^] . (5) 

(nm) 

III. HOMOGENEOUS STATES 
A. Ferromagnetic state 

In the ferromagnetic state, we have (enem) = 1 far be- 
low the Curie temperature. First, we consider the elec- 
tronic sector of the problem. The Hamiltonian Hd is 
similar to that considered in Ref. 22 and the tempera- 
ture T enters only the effective hopping integral i{T). 
To calculate the free energy of the electronic subsystem, 
we use the Hubbard I decoupling for the one-6-electron 
Green function Gb{n,n^; r - tq) = -i((rc„(r)4^^ (tq))), 
as in Ref. :22:. Here r is the time variable and T is 
the time-ordering operator. The equation of motion for 
Gf,(n, ng; r — tq) can be written in the form: 

-f^G6(n+ A,no; t - tq) + C/^(n, rig; t - rg) , (6) 

A 
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where summation is performed over sites nearest to site 
n, and Q is the 'two-particle' Green function C/(n, ng; t — 
To) = -~i{{Tcn{T)nin{T)cl^^{To))) . The equation of mo- 
tion for Q is 



d 

+ M - ) Gin, ng; r - tq) = ni5^nJ{T - tq) 



A 

+ ((f/t(r)/„+A(r)c„(r)4„(ro))) 
-((rCAW^n(r)c„(r)c1,„(ro)))} . 



(7) 



The decouphng in the first term in curly brackets gives 
((nin(r)))((f'cn+A(T)c],j,(To))) = m(Gb(n + A, no; r - 
To). Making similar decoupling in the next two terms, 
we get 

iG6(n, rig; t- tq) 
xE{(aiWWA(r))) -{{ll_^{T)Ur)))] . 

A 

In the homogeneous ferromagnetic state, the sum evi- 
dently vanishes. As a result, we obtain the closed system 
of equations for Gh and Q. In the frequency- momentum 
representation, the solutions for Gb and Q are as follows 



Gfa(k,c^) = 



CJ + /i — [/(I — Til) 
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In these expressions. 
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w = zto exp — 
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and 



c(k) = -iE' 
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where z is the number of nearest neighbors. In the case 
of the simple cubic lattice, we have z = 6 and in the 
tight-binding approximation 

C(k) = - i [cos{k^d) + cos(fc2rf) + cos(fc3d)] , (11) 

o 

where d is the lattice constant and fc* are the components 
of the wave vector. 



From Eq. ((H)), it follows that the energy spectrum of 
b electrons includes two sub-bands given by Eq. Q and 
the number of states in each sub-band depends on m . In 
the limit of large U, which is relevant to magnetic oxides, 
Eq. lO can be written as 



Ei=wC(M){l~ni), E2 = U + wCik)ni. 



(12) 



It is clear from Eq. (|12|l that the width of the lower sub- 
band is W — 2(1 — rt;)w while the width of the upper 
sub-band is 2niw. The total number of states in two sub- 
bands per site is equal to one and the number of states 
in the lower and upper sub-bands is equal to 1 — n/ and 
n;, respectively. Note that this result is valid for any 
value of U as it can be demonstrated by integration of 
the corresponding terms in the Green function Eq. 
Thus, at any doping level x and temperature T the 
upper sub-band is empty since nb + ni = n = 1 — x. In 
this case, it is reasonable to use the U ^ oo limit. The 
Green function Gb then becomes 



Gfc(k,c^) 



1 



ni 



+ /i - (1 - ni)wC{k) 



(13) 



At low temperatures, T <^ w, it is reasonable to repre- 
senr the Fermi-Dirac distribution function for b electrons 
f{E) = [exp{(£; - fi)/T} + 1]-^ by the step function 
9{fi — E). It can be shown that this approximation works 
well at t, ejT > J, and T < J, where J = zJ'S^. An 
appreciable discrepancy could arise only at doping levels 
X 1 and 1 — x <^ 1. However, at these doping lev- 
els, the homogeneous ferromagnetic state is unfavorable 
(see below) . Therefore, the number of b electrons can be 
expressed through the Green function as 



rib = -I 



duj 



d^k 



27r J (27r)= 



Gfc(k, oj + iO-sgnw)e 



(14) 



This relationship defines Ub as a function of chemical po- 
tential fi and ni. The value of ni depends on relative 
positions of /i and ejT- 

If /i < — ejT, then the Jahn- Teller-induced localization 
is unfavorable, n; = 0, rif, = 1 — x, and the chemical 
potential is found from Eq. H14|l . With the increase of 
rib, ^J■ becomes equal to — ejr, the further growth in the 
number of itinerant charge carriers is ceased, ni becomes 
nonzero, and /i is pinned at the level —ejx- In the latter 
case, the number of localized electrons can be found from 
Eq. (|14|l at /i = — ejt using the relation Ub ~ 1 — x — ni. 
As a result, we get 



1 - no 



w{l - ni) 



X, (15) 



where 



no(//) = JdE'poiE'). 



(16) 
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Figure 1: The dependence of ni (solid line) and rib (dashed 
line) on doping concentration x at ejT/w = 0.05. 



and 



C(k)) 



(17) 



is the density of states of free electrons. 

At zero doping, x = 0, the number of localized elec- 
trons ni = 1, and the bandwidth W — 0. At low doping, 
all electrons are localized, rt; = 1 — a; and W = 2xw, until 
the bottom of the band reaches the energy of I electrons 
— ejT at some critical concentration x = xi{T) = e^r/w- 
At a; > xi{T) the locahzed and band electrons coexist 
as long as x is smaller than the second critical doping 
level X2{T) at which the existence of localized electrons 
becomes unfavorable, that is, n/ = if x > X2(T). From 
Eq. we find X2{T) = 1 - 7io(-ejT/'w>). Naturally, 
such a picture can exist only for a certain relationship 
between the parameters of the model, in particular, if 
ejT < w. Since w decreases with the growth of tem- 
perature, both these critical concentrations, xi and X2, 
become larger when temperature increases. Note that 
the homogeneous ferromagnetic state can exist only at 
X > xi if Tifc ^ 0. 

The values of ti; and as the functions of doping x are 
shown in Fig.^at typical parameters of the model. These 
calculations and calculations below are performed for a 
simple cubic lattice with C(k) given by Eq. Hll|l . In this 
case, the expressions for po{E') and no(/i') in Eqs. H17() 
and 1)16(1 are found in Appendix ^ see Eqs. I|A5|1 and 
(|A6I) . respectively. 

At xi{T) < X < X2{T), the kinetic energy of h electrons 
per site can be written as 



kin 



A 



(18) 



For further purposes we introduce the function 
^(^') = 7 E<4ca) = -(1 - m)eo (-j^] , (19) 

2 . \ J. Til 



where 



dE' E'poiE') 



The kinetic energy of h electrons then reads 



Ev\ 



kin 



-wA 



(20) 



(21) 



The energy of the on-site Coulomb repulsion can be found 
using the 'two-particle' Green function Q 

(22) 

Since fi < U, the pole coming from the denominator of 
the function Q (see second Eq. ©) does not contribute 
to the integral over uj, and in the limit U — > oo, we get 



f€jT\ 

U{nwnm) ^niwA[ . (23) 

\ w / 



Note that we replace the Fermi-Dirac distribution by the 
step function. Therefore, we can omit the electron en- 
tropy term and write the free energy of electrons per site 
as a sum of Eqs. |(2U, if^ . and the JT term 



F,, = -{l-ni)wA{''-^) 
\ w / 



(24) 



For X > X2(T) when ni — 0, we should replace cjt/w by 
— /i' in formulas H21|l . (|23|l . and (|24|l . where fi' is found 
from the equation 1 — a; = no(/i'). 

Now, we consider the magnetic part of the Hamilto- 
nian, iJ,„, Eq. Following the conventional mean 

field approach for spin systems^, we replace Cm = 
{sin 0m cos 0m, sin 6'm sin 0m, cos 0m} in Eq. © by its 
mean value {0, 0, m}, where m — (cos 0m). As a result, 
-ffm decouples into the sum of N independent one-site 
Hamiltonians, = NHmo, where 



i7mo(cos0) 



-aC-^ 

V w 
-\-Jm cost 



w{T) 



(25) 



w{T) — zt{T). The value of m is determined by the 
self-consistency condition 



1 

/ 

-1 



/ due-'3«'"n(") 



(26) 



where /3 = l/T. We should also take into account that 
w is related to Hmoiu) and m by Ea. H10(l 



w{T) \ 



V2 



(27) 



/ due-'3^""(") 
-1 



6 



In order to find the Curie temperature Tc, we expand 
right hand side of Eq. (|26|l in a power series in m. Using 
Eqs. (j23l and 1(23 we find 

m ^ ai{T)l3m - a3{T)f3m^ - a^{T)[3-nt> + ... . 

The Curie temperature is found then from the condition 



ai{Tc 



w{Tc 



c, 



A 



V2" 



2V2 \ w{Tc) 



Tc ■ (28) 



If we neglect the effect of polaron band narrowing, A = 0, 
we find the exphcit expression for the Curie temperature 
Tc = [wf)A{e.jT^/ wo)/ {2^/2) - J]/3. It is clear that 
polaron band narrowing reduces Tc. 

Let us now analyze the order of the phase transition 
at T = T(7, which depends on the sign of 03 (Tc), see 
Ref. S If a3(Tc) > 0, then m tends to zero at T ^ Tc 



AT)' 



l ai{T)~T 
«3(T) 



(29) 



In this case, we have second order magnetic phase tran- 
sition. In the opposite case 03 (Tc) < 0, m behaves ap- 
proximately as {a^(T) > 0) 



m{T) 



\a3{T)\ , ai(T). 



a^iT) a3{T) 



T^Tc 



W3{Tc)\ 
asiTc) 



^0, 
(30) 

and the transition to paramagnetic (PM) state is of the 
first order. The analysis shows that a5(T) > for any 
values of parameters of the model, but the sign of a3{Tc) 
can be both positive or negative depending on the doping 
level X. At some x = X12, the coefficient 03 (Tc) changes 
its sign. At low doping, x < X12, we have a^iTc) < 
and the transition from the FM to PM state is of the 
first order. At x > X12, we have the second order phase 
transition. 

The free energy of the system per site is equal to 



^fm = Fcl — TlnZrn, Z„i — S 



AtT ^0,m^l,w 
the free energy 



Wo = zto, Z„ 



(31) 

exp(— /3J), and 



■FfmlT^o ^ -(1 - n'i)u)oA ( ^ 



In the PM state, we have m 
2S+1, and 



ejT"/ + J ■ 
0, w = w{T)/V2, Z„ 



The transition from FM to PM state does not mean 
that rif, becomes zero. As it was mentioned above, rif, ^ 
at a; > xi{T) and the value of xi increases with the tem- 
perature. At a certain temperature T*, xi{T) exceeds x. 
It is clear that T* > Tc since the FM state can exist only 
at ni, 7^ 0. Therefore, in our model, in addition to mag- 
netic phase transitions, there should exist an electronic 
transition to the state without itinerant electrons. The 
temperature T* is determined by the evident condition 
x,{T*) ^ ejTV2/w{T*) ^ X. 



B. Antiferromagnetic and canted states 

In our model, there can exist other homogeneous states 
competing with FM and PM states. It is natural to 
consider the two-sublattice antiferromagnetic and canted 
states. In the canted state, the angle v (canting an- 
gle) between the local spins belonging to two sublattices 
varies from tt (AF state) to (FM state). Here we con- 
sider AF and canting states of G type, that is, in the 
cubic lattice each site of one sublattice is surrounded by 
sites of the second sublattice. 

At finite temperatures, the canting angle is defined as 
cos = (eoGA)- In the mean field approximation used 
above, we have (eoe/^) = {eo){e/^). Thus, we can de- 
fine variable m similar to that introduced in the previous 
subsection, = cosv if j/ < 7r/2 or = — cos^ if 
f > Tr/2. In the first case, the value of m{T) and the free 
energy of the system is found in the same way as for FM 
state from Eqs. (P5|l - l|77|) . The only difference is that 
TO = 1 at T = in FM state (or = 0) while to(0) < 1 
(or v ^ 0) in the canted state. Note that the tempera- 
ture of the transition from canted to paramagnetic state, 
Tcant, is given by Eq. (|28() . where Tc should be replaced 

by Tcant- 

In the case cosj^ < 0, instead of Eq. (|25|l and H27|l . we 
have 



-1 
-1 



(32) 



/ duy/1 ~ MTOe"'^^"0(") 
wjT) -! 

w = — ^ , (33) 



V2 



where mean- field Hamiltonian Hmo{u) is given by 
Eq. H25|l . The temperature of the phase transition from 
canting to PM state is found now from the equation 



Tpm = -[l-ni] 



^T)^/ejTV2 
V2 \ w{T) ^ 



-ejTni^T\n[2S+l). 



w{Tcant) . / e,/T\/2 



2V2 



-A 



w{Tcant) 



T 



(34) 
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Figure 2: The phase diagram of the system without taking 
into account the possibility of phase separation. The regions 
1 and 2 correspond to the canted states with cos > and 
cosi/ < 0, respectively. The parameters are tjx/wo = 0.05, 
J/wo = 0.01, A = 10, and Q./wo = 0.03. 



The canted state can exist at ^ since it arises in our 
model due to the motion of conduction electrons. 

At doping levels x < xi (T), we have rtb = and the AF 
ordering is favorable. At rather low doping, x < a;i(0), 
nfe = at any temperature. In this case, the Neel temper- 
ature is independent of x and is determined by Eq. H34|) 
as = J/3. At higher doping, there can occur a tran- 
sition from the AF state to the PM state with nf, 7^ 0, 
and the Neel temperature, Tat, is determined by the com- 
parison of free energies of corresponding states. In this 
doping range, the AF state turns out to be more favor- 
able than the canted state. With the further increase 
of X, the phase with ni, has the lower energy, but 
this phase is FM rather than canted. The temperature 
of the AF-FM transition is found from the comparison 
of corresponding free energies. This transition is of the 
first order. Note that the phase diagram exhibits AF- 
FM-PM triple point a.t x — X3 {x^ < 2:12). The canted 
state can exist in the doping range 1 — x <C 1, where 
the number of itinerant electrons is too small to stabilize 
the FM state. The corresponding phase diagram in the 
X — T plane is shown in Fig. |21 We see that the compe- 
tition between the itinerant and localized electrons gives 
rise to a rather complicated magnetic phase diagram of 
the system. 



model of electronic phase separation widely discussed in 
connection with manganites and other magnetic oxides. 
Among the possible types of phase separation, we treat 
below the coexistence of different phases: AF-FM, FM- 
PM, and PM phases with different values of nt,. For sim- 
plicity, we do not include into consideration the phase 
separation involving the canted state since it exists only 
in the narrow doping range. Note that our model always 
leads to such kind of phase separation, where we have 
nfc = in one of the phases. 

We consider a system separated into two phases with 
the volume concentrations p and 1 — p. In the homo- 
geneous phases, the electron concentration per site co- 
incides with the doping level x. In the inhomogeneous 
states, the electrons can be redistributed between the re- 
gions with different phases. Let rif, 7^ in the first (F) 
phase and rif, = in the second (A) phase, the electron 
density per site in the first phase is Xf and in the second 
phase is Xa- So, doping level x lies between Xa and Xf. 
The charge conservation requires pxf + {1 — p)xa = x. 

In Fig. 121 we show the dependence of the free en- 
ergy of the most favorable homogeneous state Fhom = 
min (Ffin, Fat, ^pm, ^cant) On the doping level at different 
temperatures. TheFhom(2;) curves have two minima: one 
at a; = and another near x — X2{T). Then we could 
expect that Xa should be around zero, while Xf should 
be close to x = X2{T). 

The phase separation corresponds to the non-uniform 
charge density and we should take into account the 
Coulomb contribution to the total energy. This contri- 
bution depends on the structure of inhomogeneous state. 
To evaluate the Coulomb energy, we assume the spher- 
ical geometry of the phase-separated state. Namely, at 
p < 0.5, the sample is modelled as an aggregate of spheres 
of F phase embedded into A matrix or that of A spheres 
in F phase for p > 0.5. For this geometry, it is reasonable 
to calculate the Coulomb energy using the Wigner-Seitz 
approximation: each F or A sphere of radius Rs is sur- 
rounded by a spherical cell of radius Rceii, such as the 
volume of the cell is AttRI^h/S = V/Ng, where V is the 
volume of the sample and Ns is the number of spheres. 
The radius Rceii is related to Rg as Rg ~ P^^^Rceii for 
p < 0.5, and Rg = {1- pf'^Rceii for p > 0.5. The total 
electric charge inside this cell is zero and the Coulomb 
energy of the system is the sum of the electrostatic en- 
ergies of these cells. Following Ref. |^ we obtain the 
expression for the Coulomb energy per site Ec at p < 0.5 



IV. PHASE SEPARATION 

Up to this point, we considered only homogeneous 
states, however, it is well known that different inhomo- 
geneous states are possible in the systems with strongly 
correlated electrons. So, we should compare the free en- 
ergies of the states studied in the previous sections with 
those for inhomogeneous states. As a typical example 
of an inhomogeneous state, we analyze here the droplet 



where e is the average permittivity of the sample and d 
is the lattice constant. In the case p > 0.5, we should 
replace Xf^Xa and p ~* 1 — p. 

The b electrons in the phase-separated state are con- 
fined within a restricted volume. The corresponding size 
quantization gives rise to the change in the density of 
states. The additional contribution to the energy (per 
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Figure 3: Free energy of homogeneous state Fhom = 
min(_Ffm, -Fat, Fpm, -Fcant) VS. dopiug level x, at diflterent tem- 
peratures. The dot curve corresponds to the function X2{T). 
The parameters are ejr/wo ~ 0.05, J/wq = 0.01, A — 10, and 
n/wo = 0.03. 



site) is proportional to the total surface area between F 
and A phases, and at p < 0.5 can be written in the form 



(36) 



where surface energy (T{xf) is calculated in Appendix 1X1 
In the case p > 0.5, we should change p 1 — p. 

The Coulomb H35I) and surface H36() contributions to 
the total energy depend on the size Rs of the inhomo- 
geneities. Minimization of Ecs — Ec + with respect 
to Rs gives at p < 0.5 



Rs 



1/3 



47ru(x/ - Xa)'^ (2 - 3pi/3 +p) 



(37) 



E, 



9tt ^ 



a^a)V(x/)^ 



1/3 

p(2 



3p 



1/3 



P 



1/3 



(38) 

where u — e^/erf. 

Let us estimate the parameter u and the characteristic 
size of inhomogeneities Rg. Using typical values of pa- 
rameters for manganites z = 6, to = 0.3 eV, d = 0.4 nm, 
and e = 20 we find u « 0.18 eV, and u/wq ~ 0.15. The 
surface energy a{x) is calculated in Appendix For 
example, at x = 0.2 the optimization procedure gives 
Xa <^ Xf 0.5, p ~ 0.4 and <T{xf) w 1.4 • IO'^wq and 
from Eq. H37|l we find Rg ~ 1.5 d, that is, the inhomogene- 
ity contains Ng — 4TT{Rs/d)^ /3 « 15 unit cells. Thus, 
we see that A^^ ^ 1 for characteristic values of param- 
eters and the Wigner-Seitz approximation is applicable. 
However, this approximation overestimates the Coulomb 
contribution because of a sharp boundary of the inho- 
mogeneities. Therefore, the above values for Rg and Ng 
could be considered as lower estimates. 



In the phase-separated state, to find the values of 
and Xa at given x and T, it is necessary to minimize the 
total free energy 

FpsiXf,Xa) ^pFy{xf) + (1 - p)FA{Xa) + Ecsix f , Xa, p) 

(39) 

with respect to xj and Xa , where p = [x — Xa) / {x f — Xa) ■ 

At temperature T > T^^^ = r*(x = 1) (see Fig. 
only homogeneous PM state with rib = is possible. As 
it follows from the numerical and analytical analyzes of 
the free energy, below T^^^ , the function of two variables 
F-pa,{xf,Xa) has two minima if 2; < X2{T) for any phases 
F and A within the considered hierarchy of parameters 
(J ^ €jT < Wo)- The first minimum at the point Xa = 
Xf = X corresponds to some homogeneous state while 
the second minimum at the point Xa = 0, Xf = 2:2 (T) 
corresponds to the phase-separated state. In the limit 
u 0, the second minimum is the global minimum of the 
function Fps{xf,Xa) at < a: < 2:2 (T) and T < T^^^, 
and the PS state is favorable. When u increases, the 
range of phase separation in the plane (x, T) gradually 
narrows and disappears at some critical value Uc. Note 
that there is no localized electrons in a more metallic F 
phase (xf = X2{T)), and vice versa, in the insulating A 
phase nb{T) = since Xa = 0, as it was mentioned in 
connection to Fig. O 

Since the concentrations Xa and Xf are independent 
of the doping level, the temperatures of magnetic phase 
transitions in both phases do not depend on x. The Neel 
temperature of the AF phase is — J/3, whereas the 
Curie temperature for the FM phase can be found from 
the equation 



2V2 



e.7T 



V2' 



w{Tc) 



- J 



(40) 



The Neel and Curie temperatures in the PS state calcu- 
lated in such a way correspond to the macroscopic phases 
with the size of inhomogeneities i?s/d ^ 1. If i?s ~ rf, 
the values of Tm and Tc in the PS state can differ from 
those calculated above. 

The region where the PS state is favorable can be found 
from the comparison of free energy -Fps(a;2(7^), 0) with 
the free energies of homogeneous states at given x and 
T . The phase diagram of the system in the (x, T) plane is 
shown in Fig. 01 The range of existence for the PS state 
is bounded by the curve Ty>s{x). In the PS state, the 
content p of the metallic (rib 7^ 0) phase varies with the 
temperature and the doping level. Hence, the insulator- 
metal transition is possible when p exceeds the percola- 
tion threshold. 

Let us now discuss the transition of the system from 
the PS to a homogeneous state. The volume fraction of 
the F phase in the PS state is p{T) = x/x2{T). Depend- 
ing on the relation between the temperatures Tps, Tp, 
Tjv, and T*, the system can pass from the PS state to 
the FM [p = 1), AF (p = 0), and PM (with rif, 7^ or 
rib — 0) homogeneous states. In all cases, the number of 
itinerant electrons Ub undergoes a sudden change at the 
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Figure 4: The phase diagram of the model at ejr/wo = 
0.05, J/wo = 0.01, A = 10, n/wo = 0.03, and u/wo = 0.5. 
The numbers denote: 1. homogeneous AF phase; 2. the 
mixture of two PM states with rii, 7^ and nt = 0; 3. and 4. 
homogeneous canted states. 




Figure 5: The temperature dependence of nt &t x = 0.3, 
ejT /wo = 0.05, J/wo = 0.01, A = 10, fl/wo = 0.03, and 
u/wq = 0.5. 



transition to the homogeneous state. The temperature 
dependence of is shown in Fig. |31 



V. EFFECT OF MAGNETIC FIELD 

In this section, we consider the effect of magnetic field 
on the properties of the system. We take into account 
only the effect of the magnetic field on the local spin. 
This corresponds to the limit of classical local spin 5* 1. 
Thus, in the presence of external DC magnetic field H, we 
should add the term — MssX^n in Hamiltonian |(2Jl, 
where /is is the Bohr magneton and g is the Lande factor. 
As a result, the magnetic field term modifies only the 



magnetic Hamiltonian, Eq. 

(4cm) (t^- 



E 

(nm) 

J'5'^e„ei„ + Gnh 



(41) 



h — fiBgSU . 



In the FM state, the one-site magnetic Hamilto- 
nian H25|l corresponding to the MF approximation takes 
the form 




(42) 

where the direction of magnetic field is parallel to z axis. 
The mean value m — (Sq)/ S is found by solving the sys- 
tem of equations H26(l and H27() with Hamiltonian H42|l . 
At T < Tc, the correction to the free energy in the pres- 
ence of magnetic field SF ^ —h, whereas in paramagnetic 

phase SF h'^/Tc- 

In the AF or canted states, the result depends on 
the mutual orientation of h and the vector 1 = ((eo) — 
(ga))/2. The minimum of the free energy corresponds 
to the case h_Ll. Let the vector 1 be parallel to z axis 
and the magnetic field h be parallel to the x axis. The 
mean values of the directions of local spins in two sublat- 
tices, (gq) and (e^) can be written as (eg) = {m, 0, I}, 
(ga) = {m, 0, —I}, where m is proportional to the mag- 
netization of the system. The one-site magnetic Hamil- 
tonian then has the form 

HUO,^) ^ -A{^) (^.(T)yi±^-^) 

+Jgo(ga) — ^sin6'cos(/), (43) 
where Gq — {sin cost/), sin^sint/), cos 9} and 
eo(eA) = cos 9 + m sin 9 cos . 
The values of I and m are found from the equations 
Jdn cos0e-''^"o(«'^) 

(44) 



I 



/df^e-/5«-o(«-« 

fdfl sine cos (1)6^1^""'°'^'' 
/df2e-'3^™o(e,0) 



Expression for the effective bandwidth w (|33() now takes 
the form 



r ^r^_-fl)^,_,„rft • (45) 



V2 



At high magnetic fields h T/v ~ J, we get from 
Eq. (|44|l that ^ = 0, and the system passes from the 
AF or canted state to the FM one. The typical Neel 
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Figure 6: The temperature dependence of {nb{h) 
nb(0))/n6(0) at h/J = 0.2, x = 0.3, ejr/wo = 0.05, J/wo 
0.01, A = 10, n/wo = 0.03, and u/wo = 0.3. 



temperature in manganites is T/v ~ 100 K, and the value 
of h = J corresponds to the fields of the order of H ^ 
100 T. 

If ft-/ J <C 1, the system (I44II can be solved perturba- 
tively. The correction to the free energy of the AF state 
is 5F ^ —H'^/Tn both below and above T/v, because 
there is no spontaneous magnetization in the system at 
H = 0. The external magnetic field favors the FM state 
in comparison to the AF and canted states. In particular, 
it reduces the temperature of transition from the AF to 
FM state. The magnetic field leads to the increase in the 
effective hopping integral i due to the alignment of local 
spins and thus to the growth in the number of b electrons 
and the value of T*. At the transition from the PS to 
PM homogeneous state the magnetic field results in the 
increase of the transition temperature, and the difference 
ATps = Tps(H)-Tps(0) h. At the transition from the 
PS to FM state, ATps can be both positive and negative 
depending on the parameters. 

As was mentioned above, the number of itinerant elec- 
trons Hh differs significantly below and above Tps- There- 
fore, the shift of the transition temperature Tps with the 
magnetic field gives rise to a significant change in the 
number of itinerant electrons. The temperature depen- 
dence of the ratio {nb{H) — nb(0))/n;,(0) near the transi- 
tion temperature rps(O) is shown in Fig. El The narrow 
peak in this ratio is a manifestation of the step in nh(T) 
shown in Fig. Since the position Tps of the step in 
nh{T) depends on the magnetic field, a small change in 
H causes a significant change in the number of charge 
carriers at a fixed temperature near Tps- The number of 
itinerant charge carriers determines the value of metallic 
conductivity of the system. Thus, the large change of nt 
in magnetic field can be related to the colossal magne- 
toresistance effect. 



VI. CONCLUSIONS 

We discussed a " minimal model" dealing with the com- 
petition between the localization and metallicity in man- 
ganites. The Hamiltonian of the model takes into account 
the essential physics of strongly correlated electron sys- 
tems with the Jahn- Teller ions: it is, in fact, the Hubbard 
model with the strong electron-lattice interaction, the 
Hund's rule intraatomic coupling, and AF interatomic 
exchange between local spins. Such an approach pro- 
vides a possibility to understand the difference between 
the number of itinerant charge carriers Ub and the dop- 
ing level'^^ . It is shown that rib can be significantly lower 
than the number of the charge carriers implied by the 
doping l evel. The models of similar type were discussed 
in Refs. I26l27l3d However, the possibility of the phase 
separation was not considered in these papers. 

Here we demonstrate that in the framework of our 
model phase separation can exist in a wide range of inter- 
mediate doping concentrations disappearing at low and 
high doping level. These predictions are in agreement 
with the general features of the experimentally found 
phase diagrams of manganites^i^LS. The obtained re- 
sults suggest the existence of the droplet-type of the elec- 
tronic phase separation that was widely discussed in lit- 
erature (see, e.g., Ref.l^. We calculated the relative con- 
tent of different phases in the phase separated states and 
found the size of such droplets (ferrous). For the char- 
acteristic values of parameters, a droplet includes 10-30 
unit cells. These results could, in particular, serve as a 
key to an adequate description of the transport properties 
of manganites that could not be done in the framework 
of the single-band modelsSS;^. 

As it was mentioned above, various models related to 
the strongly correlated electrons exhibit a pronounced 
tendency toward phase separation. The effective Hamil- 
tonian of our model Q is, in fact, a generalization of the 
Falicov-Kimball modcl'^^. The latter model describes a 
system with the hybridization of an electronic band and a 
localized level. The Falicov-Kimball model is often used 
as a toy model in the analysis of heavy-fermion materials 
and it also leads to the phase separation phcnomenaii^. 
So, we believe that our approach could be applicable not 
only to manganites but also to a wider class of strongly 
correlated electron systems. Note that the analogy be- 
tween the Falicov-Kimball model and the Hamiltonian of 
the s — d type with the Jahn- Teller interaction was indi- 
cated in Ref . Isd but in the case without phase separation. 

In this paper, we analyzed the phase diagram of the 
model in the x — T plane. The effect of temperature 
manifests itself mainly in the change of effective hopping 
integral t due to the polaron band narrowing and the en- 
tropy term in the free energy due to thermal fluctuations 
of local spins. The polaron band narrowing is described 
by standard formula Q , see Refs. lslM The behavior of 
the local spins was treated using the mean field approxi- 
mation. We find that at low temperatures the system is 
in a state with a long-range magnetic order: AF, FM or 
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AF-FM phase separated state. We demonstrate that at 
high temperatures there can exist two types of the para- 
magnetic state, a usual one with nt = and that with 
Uf) 7^ 0. In the intermediate temperature range, the phase 
diagram includes different kinds of the PS states: AF- 
FM, FM-PM, and PM with different content of itinerant 
electrons. 

The applied magnetic field leads to the changes in the 
phase diagram. It evidently favors the FM ordering and, 
consequently, to the increase of the number of itinerant 
electrons. The effect of the magnetic field was analyzed 
accounting for the alignment of the local spins in the 
applied magnetic field. 

It is demonstrated that in our model the metal- 
insulator transition can take place at some characteristic 
values of the doping x corresponding to the crossover be- 
tween different kinds of the phase separation. It can be 
induced by changing temperature or the magnetic field 
and is of a percolation type. This transition can be re- 
lated to the colossal magnetoresistance effect. 

Note that in the present treatment we assume that the 
effective parameters t and ejt do not depend on the dop- 
ing level X. To verify the applicability of such approxima- 
tion, we calculated the phase diagram with t and ejT de- 
pending linearly on x. We found that the phase diagram 
remains qualitatively the same even if t and ejx vary by 
the factor 2-3 provided the hierarchy of the model pa- 
rameters ( J <C ejT < uiq) remains unchanged. 

Note also that we included to our analysis the long- 
range Coulomb interaction related to the macroscopic 
charge redistribution in the phase separated state. It 
allowed us to estimate the size of the inhomogeneities. 
At the same time, we did not take explicitly into account 
the corresponding terms in the model Hamiltonian j^J. 
However, if we would like to consider the effect of charge 
ordering, we have to include at least the nearest-neighbor 
Coulomb repulsion. 

Here we considered only the "minimal model" describ- 
ing the effect of the phase separation. Therefore, we 
did not include into consideration the possibility of the 
charge ordering to focus the discussion on the interplay 
between localized and itinerant electrons. It is not im- 
possible to include the charge ordering to the model of 
such type. The first attempt was made in Ref.|23, but the 
possibility of phase separation was not considered there. 
The effect of the charge ordering could change the results 
for X near 0.5. Therefore, it is reasonable to consider only 
the range x < 0.5, if we want to compare our predictions 
with the actual situation in the doped magnetic oxides. 
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Appendix A: SURFACE ENERGY 

In this section, we calculate the surface energy coming 
from the size quantization. The expression for the free 
energy of b and I electrons can be written in terms of the 
density of states in the form (see Eqs. l(T^ and 

^^ci = wil - ni)^ J dE' E'poiE') - ejTUi , (Al) 
-1 

where 



w{l - ni) ' 

This expression is valid for x < X2(T) and ^ = —cjt- 
The density of states for the system of itinerant electrons 
in the volume V is 

Po(ii;')-^E^(^'-C(kn)), (A2) 

n 

where momentum k varies over a discrete set of values, 
depending on boundary conditions and geometry of the 
system. The function C(kn) is normalized to unity, that 
is |C(kn)| < 1- In the thermodynamic limit V —^ oo, the 
sum in Eq. ljA2|) can be replaced by the integral over k in 
the first Brillouin zone multiplied by V/{2Tr)^. At finite 
V, we derive an approximate expression for the density of 
states in the case of cubic lattice, corresponding a small 
value of A = Sd/V, where S is the surface area and d is 
the lattice constant. 

Let the sample has the shape of a parallelepiped P 
with sides Li, L2, and L3 (in units of the lattice constant 
d). The Dirichlet boundary conditions for the conduction 
electron wave function i/'m is used, that is, Vm = for 
ni$P. In this case, the momentum kn takes values k„ = 
7rn/(La + 1), where n" = 1, 2 . . . (a = 1, 2, 3). For 
large La , we can use the trapezium rule for the sum over 

L 

E/K) = ^^Jdkf{k)-^{f{0)+f{n))+O (^-^) . 

n=l g 

As a result, for 3D sum over n we obtain 

E /(k„)^(F+|)|^/(k) (A3) 

- ^(/^[/({o,fc^fcn)+/(Kfc^fcn)] 

+/^[/({fc\o,fcn)+/({fc\^,fcn)] 
+ /^[/({fc^fc^o})+/({fe^fc^.})]) , 

where 3D and 2D momentum integrations are performed 
in the range < fc" < tt, and S = 2{LiL2 + L1L3 + L2L3) 
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is the surface area of the parallelepiped in the units of 
lattice constant. 

Using formula l|A3(l with /(k) — 5{E' — C(k)), we can 
calculate the density of states. Relation (|A3I) can be sim- 
plified for the case of cubic symmetry since C,{k^ , fc^, k^) — 



C(fc^,fc^,fc^) = C{k^ , k"^ , k^). In the absence of external 
fields, we have ((k) = C(— k), and the integration in 
Eq. jSU can be extended to — tt < fc" < tt. As a re- 
sult, we obtain for the density of states 



J 



p{E') = (i + y) poiE') - 7 / - a{oyy})) +s{E'-- canyym , (A4) 



r 



where po{E') is the density of states at 1/ — > oo. Note, 
that the density of states in the form (|A4I) depends only 
on the ratio S/V and does not depend on the shape of the 
sample. We believe, that Eq. I|A4|I is applicable for any 
geometry of the system, provided that the minimum lin- 
ear dimension L is large compared to the lattice constant 
(see, for example^Si^) . 

Let us calculate now the surface energy for the spec- 
trum in the tight-binding approximation for the simple 
cubic lattice, ((k) = — [cos(fc^d)-|-cos(fc^(i)-|-cos(A:^d)]/3. 
In the limit I^ ^ oo, the formulas for the density of 
states po{E'), density of itinerant electrons n^{p), and 
their kinetic energy eo (//) can be written in the following 
form 



the density of states Eq. (|A8I) . Note that ti; and are 
functions of A. In the considered limit A <C 1, we can 
use the perturbation technique to calculate the surface 
energy a. Representing the number of I electrons and 



the free energy F^i in the form n/ 



,(0) 



An 



(1) 



(0) 



n— S 

and HA1|I . 



Act -I- . . . 

one obtains 



and expanding Eqs. H15|) . 



\(l-nr) (AlO) 



1 - noiji'o) + fJ-'oPaipQ) 



+ C>C 



poiE') I ^j3Q)cos(i?'«), 



(A5) 





+ OC 



no 



du ( u\ sin(u) + sin (/i'lt) 



, (A6) 



du 3 (u\ 



p,' sin(/i'u) — sin(u) 



cos{p'u) — cos(u) 



where Jo is the Bessel function. Now, Eq. ljA4p can be 
rewritten as 



P{E') 



where 



1 + y )po(£^') 



A 



(A8) 



[p^;^\E' + iii)+p'i\E'-i/i)) , 



(2) 

Pi) 



+ 00 

{E')^ j ^Jl{l)cos{E'u) (A9) 



+ Ue,M-e'^o\p', + \)-^f{p',-\) (All) 



1 
3 

JO) 



4''(Mo + :t)-4'^(Mo-:t) 



(A7) where is determined by Eq. (|15|l at A = 0, /ig 



— ejT/['w(l — n|°^)], n'^\p') and £^^\p') are the number 
of b electrons and their energy in 2D case at A = 0, 



+ 00 



n^^\p') 



du 2 (u\ 
~ V 3 J 



u\ sm I 



+ sin{p'u) 



(A12) 



p' sm{p'u) 



COs{p'u) — COS (^) 



(A13) 



is the density of states in the 2D case. The number n; and 
the free energy of electrons (in the case n/ ^ 0) is given by 
Eqs. H15|l and (jAl|) . where instead of po{E') we should use 



At doping concentration x > X2{T) when n; = 0, we 
should use the equation 1 — x — no{p') for the chemical 
potential, where Pq{E') is substituted by p{E'). As a 
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Figure 7: The surface energy vs doping concentration x at 
ejT/w = 0.1. 



result, we obtain the expression for the surface energy 



^ n^o\f^o + ^)+ n^Pit^'o -^)- ^Mf^'o) 



£o(/^o) 



1 



(AM) 



where fi^ is found from the equation 1 ~ x — no(/io) 
at A = 0. The dependence a{x) is shown in Fig. |7| 
The function cr{x) is discontinuous at a; = X2{T). This 
singularity stems from the kink in the free energy Fci{x) 
(see Fig. El. 

In the approximation under discussion, the corre- 
sponding corrections to the magnetic contribution to the 
free energy Fm are of the order of A^, and therefore we 
omit them. 
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